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Abstract 

We study the photon counting noise in optical interferometers used for gravitational wave detec- 
tion. In order to reduce quantum noise a squeezed vacuum is injected into the usually unused input 
port. It is investigated under which conditions the gravitational wave signal may be amplified with- 
out increasing counting noise concurrently. Such a possibility was suggested as a consequence of the 
entanglement of the two output ports of a beam splitter. We find that amplification without con- 
current increase of noise is not possible for reasonable squeezing parameters. Photon distributions 
for various beam splitter angles and squeezing parameters are calculated. 

PACS numbers: 42.50.Ar, 42.50.Ex, 42.50. St 
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I. INTRODUCTION 



Optical interferometers for the detection of gravitational waves need an extremely accu- 
rate control of various sources of noise. In an effort to reduce quantum-mechanical noise 
in such interferometers, Caves proposed the squeezed state technique (!']: into the normally 
unused port of the interferometer a squeezed vacuum state is injected. Details of this tech- 
nique have since been analyzed in many experimental and theoretical investigations (e.g. 
Refs. |2, 3|). 

A central element of an interferometer is a beam splitter, and as can be easily shown, a 
photon state, which is a product of the states in each input port, is not a product of the 
photon states in each output port of the beam splitter. Instead, the two output states of 
a beam splitter are entangled, and it was the subject of recent work by Barak and Ben- 
Aryeh (4] to investigate the consequences of this entanglement for the photon statistics of 
an optical interferometer. 

In particular, it was suggested that under certain conditions the gravitational wave signal 
may by amplified significantly without a corresponding increase in counting noise. This effect 
was attributed to the entanglement effects mentioned above in connection with squeezing. 

It is the purpose of this paper to investigate this surprising proposal in detail. To this end 
we determine the photon distributions in the output state of a beam splitter for weak and 
strong coherent states injected into one of the input ports of the beam splitter in addition 
to the squeezed vacuum in the other. Our results are in disagreement with Ref. {4] for both 
weak and strong input states. In particular, we cannot confirm the finding in Ref. (4] that 
the gravitational wave signal may be amplified without concurrent increase of noise. 

In order to handle the entanglement effects one needs to disentangle exponential oper- 
ators. This is rather simple in the present case, and we describe a Lie algebraic method, 
where the disentangling coefficients are calculated numerically from a set of coupled nonlin- 
ear equations. We provide a complete graphical overview of the disentangling coefficients 
as a function of the squeezing factor of the injected squeezed vacuum state and the angle of 
the beam splitter. 

In section 2 we develop a general formula for the calculation of the photon number 
distributions in the two output ports of a beam splitter. Furthermore, we discuss the special 
case where the beam splitter is oriented close to 90° to an incoming strong coherent state 
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with a squeezed vacuum entering the other port. Numerical results will be discussed in 
section 3. A brief summary concludes the paper. Technical details are relegated to a few 
Appendices. 



II. PHOTON STATISTICS IN A MICHELSON INTERFEROMETER 

We consider a beam splitter and inject a coherent state into port 1 and a squeezed vacuum 
state into port 2. The incoming photon state is therefore described by 

foMa,C)> = S 2 (()Di («)|0,0) (1) 



with 

Di(a) = exp (aa\ - a*ai) , S 2 (() = exp ( y d 2 , - ^dfj ■ ( 2 ) 

The coherence parameter a and the squeezing parameter ( are complex numbers. The 
photon field operators dj and d] fulfill the boson commutation relation [a*, dj] = 5ij. 

After passing the beam splitter the field is described by the rotated field operators h and 

b\ with [k, St] = (j w Q 

/ 



Oi 
«2 
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(3) 



cos 7 sin 7 
y — sin 7 cos 7 

The parameter 7 parameterizes the splitting ratio of the beam splitter with respect to the 
incoming beam in port 1. 

Now expressing the a operators in terms of the b operators enables us to write the photon 
state after passing the beam splitter as follows 

|Vw(«,C,7)> =exp(|C|A) J D 1 (acos7) J D 2 (asin7)|0,0) (4) 

with 

A = §i sin 2 7 + §2 cos 2 7 + §12 sin 7 cos 7 (5) 

and 

§t = 2e (c *^ 2 " c ^ t2) Su = ^ (c ^ 2 " Chb2) - (6) 

In the absence of the operator Su the output would be obviously characterized by squeezed 
coherent states in both output ports. However, these states are entangled via §12 and the 
output state cannot be factorized. This fact significantly complicates evaluation of the 
photon statistics of the output state. 



However, it is possible to use Lie algebraic disentangling techniques in order to rewrite 
the output state in a way which enables the determination of photon distributions. To this 
end we first define the additional operator i 12 , 

ii2 = hbi - b\b 2 . (7) 

The operators sj, §2, §12, and £ 12 form a closed Lie algebra with the commutation relations, 

[si, s 2 ) = 0, [sia, i 12 ] = -2Si + 2s 2 , 

[Sl, S12] = *12, [S2, S12] = -iu, (8) 
[Sl,il2] — S12, [52,^12] — ~~ S12. 

As a consequence it is possible to write Eq. ([9} as follows js, 7] 

l^out) = exp(a T ^2)exp(asSi2)exp(aiSi)exp(CT2S 2 )-Di(acos7)£)2(asin7)|0,0). (9) 

The output state is now expressed in terms of two squeezed coherent states entangled via 
the operators exp(a"Tti2) and exp(<7s,§i 2 ). The coefficients <jt, &s, ffi, &2 are real functions 
of the input parameters r = |£| and 7. A simple method for the numerical determination of 
these parameters is described in Appendix [A] 

In order to determine the photon statistics of the output state we need to determine its 
number (Fock) representation. We start out from the number representation of a squeezed 



coherent state js] 



^(CW)|o> = x;-rr/»(C,/3)|n> (10) 

n=o Vnl 

with 



x os (e i0 tanhrW 2 / 1,.„ |2 l0 n2 , \ / /3e~* e/2 \ 

£ = re* 9 , and if n the Hermite polynomials. Now applying the two exponential operators 
exp(cry£i 2 ) and exp (a sS 12) on a product of such states yields the Fock representation of 
the output state. After some straightforward but tedious algebra one obtains the photon 
distribution P ni „ 2 in the two output ports of the interferometer 
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(12) 



with 



Ni = rix - m 2 - m 3 + 7714, N 2 = n 2 - m 2 + m 3 - m.4, 

M 2 = min(ni — m 3 + m 4 , n 2 + m 3 — m 4 ), (13) 

A ™ m3m4 = ^A^Vt^t 4 (ni + m 4 )! (n 2 + m 3 - m 4 )! x 
ni ™ 2 mi!m2!r?T.3!m 4 ! (n 2 — m 4 )! jVilA^! 

e /is(l+ni+n2-2m2)/2 e ^ T (-«i+n2-2m4)/2 

( = re ld , /3i = a;cos7, /32 = asin7, 

A5 = e* e tanh o"5, ^5 = — 2 log(cosh erg)), ^5 = e~ %e tanh. as, 
\ T = tan ax, Ht = — 21og(coso"T)), = tancr T . 

The sum over mi in Eq. f lT2|) is in principle unrestricted above. In practice a suitable upper 
limit must be chosen such that the probability is correctly normalized. Most of the numerical 
results in the following are calculated from Eq. ( TT2l . A few details used for the derivation 
of Eq. (112j) are given in Appendix [B] 

In addition to the general case discussed above, we investigate the special case 7 = tt/2 + 5 
for which one finds (see Appendix lAl <j\ = r, o 2 = 0, as = — 5sinhr, and ax = 5(1 — cosh r), 
if 5 is sufficiently small. We furthermore assume a very strong coherent state incoming in 
port 1, such that the b 2 and b 2 operators can be replaced in the entanglement factors in 
Eq. © by their expectation values a and a*, respectively. The output state can then be 
written as 

|Vw) = Di(-a5(l -coshr))A(-5a*e ie sinhr)S 1 (C) J Di(-«5) J D 2 (a)|0,0) (14) 



with £ = re . A similar result was obtained along somewhat different lines in Ref. 
however with opposite sign in the first factor. The operators with index 1 may be combined 
using the relations given in Appendix [C] and we obtain for the output state 

IVW) = e^^StiOH-^D^MO) (15) 

with a = \a\e l< ^ and 

n(r,9-2<p) = \K\e iX = coshr + e i(e_20) sinh r, (16) 
A(r,6,cf>) = -^sin(fl-20)e- 2 ^sinh(2r). 

As one can see from Eq. (fl5|) . a strong coherent state with coherence parameter a exits 
through port 2 of the interferometer and a weak squeezed coherent state with coherence 



parameter — adn and squeezing parameter ( exits through port 1. Note, that the coherence 
parameter depends on the squeezing parameter. The phase factor e 4 M 2(52A i s irrelevant for 
the determination of the photon statistics. 

The photon statistics in port 1 is immediately obtained from Eq. f[T0"j) 

P ni = -^\f ni (C-a5K)\ 2 . (17) 
nil 

The mean and the variance of this distribution may be obtained analytically js] 

(m) = <5 2 |a| 2 (cosh2r + cos(# - 2(f)) sinh2r) x 

(cosh2r — cos(# — 20 — 2A)) sinh2r) + sinh 2 r, 
(An 2 ) = <5 2 |a| 2 (cosh2r + cos(#-20)sinh2r) x (18) 
(cosh4r — cos(6 l — 20 — 2A)) sinh4r) + 2 sinh 2 r cosh 2 r. 

The phase A(r, 9 — 20) of k defined in Eq. ( fl6l) depends on the squeezing factor r and the 
angle 9 — 20. This dependence complicates the interpretation of these formulas. We will 
discuss results calculated from Eqs. f lT7|) and ffTB"]) in section IIIIBI 

III. NUMERICAL RESULTS 

As is obvious from Eqs. (Q and (Tl2|) . in order to determine the photon statistics of the 
output state one needs to calculate the disentangling coefficients <Ji, cr 2 , &s, and ax- They 
are easily obtained as solutions of the nonlinear equations (1A6j) derived in Appendix lAl 
Numerical results are plotted in Fig. [TJ The figures clearly show the special values <j\ — r, 
02 = , as = , and <7t = for 7 = ir/2, and o\ = r/2, cr 2 = r/2, as = r/2, and o T = 
for 7 = 7r/4. We consider squeezing factors up to r = 1.5 in our numerical work, since 
at the present time the largest squeezing factor experimentally realized is about r = 1.3 
corresponding to a maximum squeezing of about -11.5 dB J9]. 

A. Weak coherent field and 7 = it /A and 7 = it/ 8 

In Fig. [2] we plot photon distributions in the output port 1 of a beam splitter calculated 
from Eq. (Tl2|) for a weak coherent field (|a| 2 = 20) injected into input port 1. The beam 
splitter is set at an angle 7 = 7r/4, and squeezed vacua with various squeezing parameters r 
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FIG. 1: (Color online) Disentangling coefficients a±, 02, as, and gt as a function of the beam 
splitter angle 7 and squeezing factor r. 

are injected into port 2. For 7 = ir/4 one finds analytically o\ = r/2, o~ 2 = r /2, erg = r/2, 
and cry = 0, and as a consequence, the entanglement operator e artl2 equals unity. We observe 
that squeezing slightly reduces the variance of the distributions in the output port up to 
about r = 0.7, and that it slightly shifts the mean of the distributions to larger photon 
numbers. This particular case was studied before by Barak and Ben-Aryeh 4|, and we note 
that our results are in contrast to Ref. J<4] where a down shift of the mean was predicted. (We 
also cannot confirm that the probabilities calculated in Ref. 4] are correctly normalized.) 
We checked, that the reasons for the discrepancies between our results and Ref. |4j are 
missing faculty symbols in all factors under the square roots in equations (36) and (39) of 



Ref. 



4]. They arise if one applies a power of a creation or destruction operator on a Fock 
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FIG. 2: (Color online) Photon distributions for \a\ 2 = 20, 7 = vr/4, and various squeezing factors 
r = 0, 0.3, 0.6, 0.9, 1.2, 1.5. On the right of the figure we plot the mean and the variance squared 
of the distributions as a function of the squeezing factor r (full lines). The dashed lines show 
calculations without the entangling factors. 

state. 

The results shown Fig. [2] clearly indicate that in order to determine the photon statistics 
of the output of a beam splitter the entanglement factors must not be neglected. In fact, 
calculations without the entanglement factors would predict an improvement of the counting 
statistics by injecting a squeezed vacuum, which is not justified. Nevertheless there is an 
optimum squeezing parameter of about r = 0.7, where we find a small improvement of the 
counting statistics by squeezing. 

Fig. US shows similar results for a beam splitter set at an angle 7 = ir/8 and otherwise the 
same parameters. For this angle we need the full numerical solution for the disentangling 
coefficients as discussed above, in particular, in this case both entangling operators e artl2 and 
e v.ssi2 contribute. Here, we see squeezing essentially has no effect on the distributions. Both 
the mean and the variance cannot be significantly influenced by squeezing. Figure [3] shows 
that calculations without the entanglement factors produce a completely wrong indication 
of an improved photon statistics gained by injecting a squeezed vacuum. 
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FIG. 3: (Color online) Photon distributions for \a\ 2 = 20, 7 = vr/8, and various squeezing factors 
r = 0, 0.3, 0.6, 0.9, 1.2, 1.5. On the right of the figure we plot the mean and the variance squared 
of the distributions as a function of the squeezing factor r (full lines). The dashed lines show 
calculations without the entangling factors. 

B. Strong coherent field and 7 = tt/2 ± 5 



In our second numerical example, we study a strong coherent field injected into port 1 
of the beam splitter. According to Eq. (1141) the output in port 2 is then a strong coherent 
state, and the output in port 1 (dark port) a weak squeezed coherent state. Its coherence 
parameter — aSn(Q depends on the squeezing parameter ( of the squeezed vacuum state 
injected into input port 2. In the Fig. H]the dependence of k on the squeezing parameter ( 
is shown. A significantly stronger increase of k with increasing r was predicted in Ref. jj]. 
The difference between the results presented here and those obtained in Ref. {4] stem from 
the sign error pointed out after Eq. (|T4l) . 

With k as input we calculate the photon number distribution for the output state from 
Eqs. (|T7j) for \5a\ 2 = 20 and \5a\ 2 = 500, respectively. Furthermore, we determine the mean 
as well as the variance of these distributions from Eq. ( fl~8l) . Results of these calculations are 
shown in Figs. [5] and [6j 

Optimally we choose a squeezing parameter ( and a coherence parameter a so as to 
amplify the signal in the output port 1 without increasing its noise. Unfortunately, this 
does not seem 10 oe p OSSlb ,e _ to the c 0nclusl0ns reach ed in Ref. Q. jjy* of 
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FIG. 4: (Color online) Modulus \k\ and phase A of the coherence parameter k = \n\e = coshr + 
e i(0-24>) gjj^j- as a function of the squeezing factor r. The label of the curves indicates the angle 
6 - 2cp. 

0.20 1 ■ . . ■ , ■ . . . , ■ . . . , ■ ■ . . 1 




FIG. 5: (Color online) Photon number probability distribution for |a<5| 2 = 20, 6 = 2(f), and different 
values of the squeezing parameter r = 0, 0.3, 0.6, 0.9, 1.2, 1.5 in the dark port. On the right 
hand side of the plot we show the mean and variance squared of these distributions. 

Eqs. ( 1T8|) as a function of r and 9 — 2(f) indicates that a significant amplification of the signal 
in output port 1 is not possible within the parameter ranges experimentally accessible and 
considered here. In particular, for 9 = 2(f) one finds A = and (ni) = <5 2 1 ct | 2 + sinh 2 r. 

Only the noise can be slightly influenced: a numerical study of the variance of the pho- 
ton distribution given in Eqs. ([IE]) shows that 9 = 20 and A = is the optimal choice. 
Then, one finds (An 2 ) = <5 2 |a| 2 e _2r + 2 sinh 2 r cosh 2 r. The minimum of the variance de- 
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FIG. 6: (Color online) Photon number probability distribution for |a<5| 2 = 500, 6 = 2(j>, and 
different values of the squeezing parameter r = 0, 0.3, 0.6, 0.9, 1.2, 1.5 in the dark port. On the 
right hand side of the plot we show the mean and the variance squared of these distributions. 
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FIG. 7: (Color online) Optimal squeezing parameter as a function of |a<5| 2 determined from the 
minimum of the variance of the distribution function. 

pends on the parameter ad, and in Fig. [7] we show the optimal squeezing parameter as a 
function of \5a\ 2 . This squeezing parameter should be chosen in order to minimize photon 
counting uncertainties. For large squeezing parameters the distributions show characteristic 
oscillations 



IV. SUMMARY AND CONCLUSIONS 



In a recent paper by Barak and Ben-Aryeh 4j it was suggested that due to entanglement 
effects in the output state of a beam splitter, the gravitational wave signal observed with 
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an interferometer may be amplified without simultaneously increasing quantum mechanical 
photon counting noise. 

In this paper we reanalyze such entanglement effects, and cannot confirm the significant 
signal amplification predicted in Ref. ^| but agree that consideration of such entanglement 
effects is important in order to calculate the photon counting noise in optical interferometers. 
We observe that for a strong coherent input in port 1 of the beam splitter and a squeezed 
vacuum injected in port 2 photon counting statistics slightly improves if the beam splitter 
is oriented slight off the right angle with the incoming strong coherent field. The differences 
between the calculations presented here and those of Ref. J4] can be traced to algebraic 
issues, which are discussed explicitly in section IIHI 

It appears that entanglement effects must be carefully studied for any particular interfer- 
ometer design, a task for which the formulas developed in the present paper may be helpful, 
in particular the general result Eq. ( !T2|) for the calculation of photon distributions in the 
output state. 



Appendix A: Disentangling 



In order to solve the disentangling problem (r = 

exp(rA) = exp(a r ti 2 ) exp(<T5Si 2 ) exp(aiSi) exp(cr 2 s 2 ) 



(Al) 



with A given in Eq. ([5]), we consider the Lie algebra of the operators jsi, s 2 , ,§i 2 , ti 2 | defined 
in Eqs. and ([7]) with their commutators given in Eq. ([8]). The corresponding matrix 
representations of these operators 

/ nn \ / 



si 




S-2 




812 




t 
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fulfill the same commutation relations. 

The matrix representation of the operator A in Eq. f lAlj) is then obtained as 



.4 



^ sin 2 7 



cos 7 sin 7 



l — cos 7 sin 7 cos 7 

Since A n = A for all n G N we find for the left hand side of Eq. (IA1I) 

(rA) n 



exp(ryl) = ^2 



n=0 



I + (e r -l)A. 



(A2) 



(A3) 



(A4) 
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The matrix representation of the right hand side Eq. (lAlj) is also easily calculated 

exp(CT T Ti2)exp(cr 5 si2)exp(aiSi)exp(a 2 ,S2) = (A5) 
' e ai (sin o"t sinh erg + cos Of cosh as) — e " 2 (sin <jt cosh as + cos Of sinh as) 
. e"" 1 (sin <tt cosh (T5 — cos ot sinh 05) — e°" 2 (sin ot sinh 05 — cos ar cosh cs) 

By equating Eqs. ( 1A4j) and flA5[) one obtains four equations for the four parameters 
o"t, ers, 01, o"2 to be determined 

e r sin 2 7 + cos 2 7 = e' 71 (sin <tt sinh as + cos o"t cosh as), 

(1 — e r ) COS7 sin 7 = — e°" 2 (sin o"t cosh cr 5 + cos <jt sinh cr^), (A6) 

(1 — e r ) cos 7 sin 7 = e CTl (sin ot cosh as — cos or sinh as) , 

e r cos 2 7 + sin 2 7 = — e CT2 (sin <tt sinh as — cos a? cosh 0-5). 

This set of equations is easily solved numerically. Results are shown in Fig. [TJ 
For 7 = Ti / 2 solutions of Eqs. (IA6I) can easily be found analytically, 

7 = tt/2, a x = r, a 2 = 0, a 5 = 0, a T = 0. (A7) 

If one expands the trigonometric functions on the right hand side of Eqs. (IA6I) keeping only 
linear terms in the disentangling coefficients, one finds 

7 = 7r/2 + 5, oi = r, a 2 = 0, as = — <5sinhr, ot = 5(1 — cosh r). (A8) 

For 7 = 7r/4 one finds 

7 = vr/4, en = r/2, a 2 = r/2, cr s = r/2, a T = 0. (A9) 

Appendix B: Evaluation of the entanglement factors 

The entanglement factors exp^axiu) and exp(crg,Si2) must be further transformed for 
convenient calculations. Using the fact that the operators b\b\, 6162, and |(1 + b\bi + 62^2) 
form an SU(1,1) Lie algebra one easily finds that 

exp[a"5Si2] = exp (bible 10 tanh as^j exp (— (1 + hi + h 2 ) ln(cosh as)) exp (— &i& 2 e -ie tanh as^j ■ 

(Bl) 
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as was proved e.g. in Ref. [lo| . Similarly one proves that 

exp[aT^i2] = exp (bible 18 tan a^j exp (— (ni — n 2 ) ln(cos cr T )) exp fej& 2 e~ t6 ' tan a^j • (B2) 

using the fact that the operators bibl, b\b 2 , and — | (6161 — 62^2) form an SU(2) Lie algebra. 
Expanding out the exponentials in the expressions above allows for convenient application 
of these operators on Fock states. 

Appendix C: Transformations 

The following transformations for products of displacement operators and squeezing op- 
erators can be easily proved from the Baker-Campbell-Haussdorff theorem and the Boson 
commutation relations, 



D(a2)D(ai) = D(a>\ + a 2 ) exp 

D(a)S(C) = S(()D(a coshr + a*e*sinhr), ( = re ie . (CI) 



These relations are used to obtain Eq. ( 115 
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